Mono-allelic KCNB2 variants lead to a neurodevelopmental syndrome caused by altered channel inactivation.

Ion channels mediate voltage fluxes or action potentials that are central to the functioning of excitable cells such as neurons. The KCNB family of voltage-gated potassium channels (Kv) consists of two members (KCNB1 and KCNB2) encoded by KCNB1 and KCNB2, respectively. These channels are major contributors to delayed rectifier potassium currents arising from the neuronal soma which modulate overall excitability of neurons. In this study, we identified several mono-allelic pathogenic missense variants in KCNB2, in individuals with a neurodevelopmental syndrome with epilepsy and autism in some individuals. Recurrent dysmorphisms included a broad forehead, synophrys, and digital anomalies. Additionally, we selected three variants where genetic transmission has not been assessed, from two epilepsy studies, for inclusion in our experiments. We characterized channel properties of these variants by expressing them in oocytes of Xenopus laevis and conducting cut-open oocyte voltage clamp electrophysiology. Our datasets indicate no significant change in absolute conductance and conductance-voltage relationships of most disease variants as compared to wild type (WT), when expressed either alone or co-expressed with WT-KCNB2. However, variants c.1141A>G (p.Thr381Ala) and c.641C>T (p.Thr214Met) show complete abrogation of currents when expressed alone with the former exhibiting a left shift in activation midpoint when expressed alone or with WT-KCNB2. The variants we studied, nevertheless, show collective features of increased inactivation shifted to hyperpolarized potentials. We suggest that the effects of the variants on channel inactivation result in hyper-excitability of neurons, which contributes to disease manifestations.

TREX tetramer disruption alters RNA processing necessary for corticogenesis in THOC6 Intellectual Disability Syndrome.

THOC6 variants are the genetic basis of autosomal recessive THOC6 Intellectual Disability Syndrome (TIDS). THOC6 is critical for mammalian Transcription Export complex (TREX) tetramer formation, which is composed of four six-subunit THO monomers. The TREX tetramer facilitates mammalian RNA processing, in addition to the nuclear mRNA export functions of the TREX dimer conserved through yeast. Human and mouse TIDS model systems revealed novel THOC6-dependent, species-specific TREX tetramer functions. Germline biallelic Thoc6 loss-of-function (LOF) variants result in mouse embryonic lethality. Biallelic THOC6 LOF variants reduce the binding affinity of ALYREF to THOC5 without affecting the protein expression of TREX members, implicating impaired TREX tetramer formation. Defects in RNA nuclear export functions were not detected in biallelic THOC6 LOF human neural cells. Instead, mis-splicing was detected in human and mouse neural tissue, revealing novel THOC6-mediated TREX coordination of mRNA processing. We demonstrate that THOC6 is required for key signaling pathways known to regulate the transition from proliferative to neurogenic divisions during human corticogenesis. Together, these findings implicate altered RNA processing in the developmental biology of TIDS neuropathology.

Variants in the WDR44 WD40-repeat domain cause a spectrum of ciliopathy by impairing ciliogenesis initiation.

WDR44 prevents ciliogenesis initiation by regulating RAB11-dependent vesicle trafficking. Here, we describe male patients with missense and nonsense variants within the WD40 repeats (WDR) of WDR44, an X-linked gene product, who display ciliopathy-related developmental phenotypes that we can model in zebrafish. The patient phenotypic spectrum includes developmental delay/intellectual disability, hypotonia, distinct craniofacial features and variable presence of brain, renal, cardiac and musculoskeletal abnormalities. We demonstrate that WDR44 variants associated with more severe disease impair ciliogenesis initiation and ciliary signaling. Because WDR44 negatively regulates ciliogenesis, it was surprising that pathogenic missense variants showed reduced abundance, which we link to misfolding of WDR autonomous repeats and degradation by the proteasome. We discover that disease severity correlates with increased RAB11 binding, which we propose drives ciliogenesis initiation dysregulation. Finally, we discover interdomain interactions between the WDR and NH2-terminal region that contains the RAB11 binding domain (RBD) and show patient variants disrupt this association. This study provides new insights into WDR44 WDR structure and characterizes a new syndrome that could result from impaired ciliogenesis.

Diagnostic utility and reporting recommendations for clinical DNA methylation episignature testing in genetically undiagnosed rare diseases.

PURPOSE: This study aims to assess the diagnostic utility and provide reporting recommendations for clinical DNA methylation episignature testing based on the cohort of patients tested through the EpiSign Clinical Testing Network. METHODS: The EpiSign assay utilized unsupervised clustering techniques and a support vector machine-based classification algorithm to compare each patient's genome-wide DNA methylation profile with the EpiSign Knowledge Database, yielding the result that was reported. An international working group, representing distinct EpiSign Clinical Testing Network health jurisdictions, collaborated to establish recommendations for interpretation and reporting of episignature testing. RESULTS: Among 2399 cases analyzed, 1667 cases underwent a comprehensive screen of validated episignatures, imprinting, and promoter regions, resulting in 18.7% (312/1667) positive reports. The remaining 732 referrals underwent targeted episignature analysis for assessment of sequence or copy-number variants (CNVs) of uncertain significance or for assessment of clinical diagnoses without confirmed molecular findings, and 32.4% (237/732) were positive. Cases with detailed clinical information were highlighted to describe various utility scenarios for episignature testing. CONCLUSION: Clinical DNA methylation testing including episignatures, imprinting, and promoter analysis provided by an integrated network of clinical laboratories enables test standardization and demonstrates significant diagnostic yield and clinical utility beyond DNA sequence analysis in rare diseases.

Spliceosome malfunction causes neurodevelopmental disorders with overlapping features.

Pre-mRNA splicing is a highly coordinated process. While its dysregulation has been linked to neurological deficits, our understanding of the underlying molecular and cellular mechanisms remains limited. We implicated pathogenic variants in U2AF2 and PRPF19, encoding spliceosome subunits in neurodevelopmental disorders (NDDs), by identifying 46 unrelated individuals with 23 de novo U2AF2 missense variants (including 7 recurrent variants in 30 individuals) and 6 individuals with de novo PRPF19 variants. Eight U2AF2 variants dysregulated splicing of a model substrate. Neuritogenesis was reduced in human neurons differentiated from human pluripotent stem cells carrying two U2AF2 hyper-recurrent variants. Neural loss of function (LoF) of the Drosophila orthologs U2af50 and Prp19 led to lethality, abnormal mushroom body (MB) patterning, and social deficits, which were differentially rescued by wild-type and mutant U2AF2 or PRPF19. Transcriptome profiling revealed splicing substrates or effectors (including Rbfox1, a third splicing factor), which rescued MB defects in U2af50-deficient flies. Upon reanalysis of negative clinical exomes followed by data sharing, we further identified 6 patients with NDD who carried RBFOX1 missense variants which, by in vitro testing, showed LoF. Our study implicates 3 splicing factors as NDD-causative genes and establishes a genetic network with hierarchy underlying human brain development and function.

Further delineation of the rare GDACCF (global developmental delay, absent or hypoplastic corpus callosum, dysmorphic facies syndrome): genotype and phenotype of 22 patients with ZNF148 mutations.

BACKGROUND: Pathogenic variants in the zinc finger protein coding genes are rare causes of intellectual disability and congenital malformations. Mutations in the ZNF148 gene causing GDACCF syndrome (global developmental delay, absent or hypoplastic corpus callosum, dysmorphic facies; MIM #617260) have been reported in five individuals so far. METHODS: As a result of an international collaboration using GeneMatcher Phenome Central Repository and personal communications, here we describe the clinical and molecular genetic characteristics of 22 previously unreported individuals. RESULTS: The core clinical phenotype is characterised by developmental delay particularly in the domain of speech development, postnatal growth retardation, microcephaly and facial dysmorphism. Corpus callosum abnormalities appear less frequently than suggested by previous observations. The identified mutations concerned nonsense or frameshift variants that were mainly located in the last exon of the ZNF148 gene. Heterozygous deletion including the entire ZNF148 gene was found in only one case. Most mutations occurred de novo, but were inherited from an affected parent in two families. CONCLUSION: The GDACCF syndrome is clinically diverse, and a genotype-first approach, that is, exome sequencing is recommended for establishing a genetic diagnosis rather than a phenotype-first approach. However, the syndrome may be suspected based on some recurrent, recognisable features. Corpus callosum anomalies were not as constant as previously suggested, we therefore recommend to replace the term 'GDACCF syndrome' with 'ZNF148-related neurodevelopmental disorder'.

Recurrent de-novo gain-of-function mutation in SPTLC2 confirms dysregulated sphingolipid production to cause juvenile amyotrophic lateral sclerosis.

BACKGROUND: Amyotrophic lateral sclerosis (ALS) leads to paralysis and death by progressive degeneration of motor neurons. Recently, specific gain-of-function mutations in SPTLC1 were identified in patients with juvenile form of ALS. SPTLC2 encodes the second catalytic subunit of the serine-palmitoyltransferase (SPT) complex. METHODS: We used the GENESIS platform to screen 700 ALS whole-genome and whole-exome data sets for variants in SPTLC2. The de-novo status was confirmed by Sanger sequencing. Sphingolipidomics was performed using liquid chromatography and high-resolution mass spectrometry. RESULTS: Two unrelated patients presented with early-onset progressive proximal and distal muscle weakness, oral fasciculations, and pyramidal signs. Both patients carried the novel de-novo SPTLC2 mutation, c.203T>G, p.Met68Arg. This variant lies within a single short transmembrane domain of SPTLC2, suggesting that the mutation renders the SPT complex irresponsive to regulation through ORMDL3. Confirming this hypothesis, ceramide and complex sphingolipid levels were significantly increased in patient plasma. Accordingly, excessive sphingolipid production was shown in mutant-expressing human embryonic kindney (HEK) cells. CONCLUSIONS: Specific gain-of-function mutations in both core subunits affect the homoeostatic control of SPT. SPTLC2 represents a new Mendelian ALS gene, highlighting a key role of dysregulated sphingolipid synthesis in the pathogenesis of juvenile ALS. Given the direct interaction of SPTLC1 and SPTLC2, this knowledge might open new therapeutic avenues for motor neuron diseases.

Chimeric Antigen Receptor T Cell Therapy Targeting Epithelial Cell Adhesion Molecule in Gastric Cancer: Mechanisms of Tumor Resistance.

Epithelial cell adhesion molecule (EpCAM) is a tumor-associated antigen that is frequently overexpressed in various carcinomas. We have developed chimeric antigen receptor (CAR) T cells specifically targeting EpCAM for the treatment of gastric cancer. This study sought to unravel the precise mechanisms by which tumors evade immune surveillance and develop resistance to CAR T cell therapy. Through a combination of whole-body CAR T cell imaging and single-cell multiomic analyses, we uncovered intricate interactions between tumors and tumor-infiltrating lymphocytes (TILs). In a gastric cancer model, tumor-infiltrating CD8 T cells exhibited both cytotoxic and exhausted phenotypes, while CD4 T cells were mainly regulatory T cells. A T cell receptor (TCR) clonal analysis provided evidence of CAR T cell proliferation and clonal expansion within resistant tumors, which was substantiated by whole-body CAR T cell imaging. Furthermore, single-cell transcriptomics showed that tumor cells in mice with refractory or relapsing outcomes were enriched for genes involved in major histocompatibility complex (MHC) and antigen presentation pathways, interferon-γ and interferon-α responses, mitochondrial activities, and a set of genes (e.g., CD74, IDO1, IFI27) linked to tumor progression and unfavorable disease prognoses. This research highlights an approach that combines imaging and multiomic methodologies to concurrently characterize the evolution of tumors and the differentiation of CAR T cells.

CAR+ and CAR- T cells share a differentiation trajectory into an NK-like subset after CD19 CAR T cell infusion in patients with B cell malignancies.

Chimeric antigen receptor (CAR) T cell therapy is effective in treating B cell malignancies, but factors influencing the persistence of functional CAR+ T cells, such as product composition, patients' lymphodepletion, and immune reconstitution, are not well understood. To shed light on this issue, here we conduct a single-cell multi-omics analysis of transcriptional, clonal, and phenotypic profiles from pre- to 1-month post-infusion of CAR+ and CAR- T cells from patients from a CARTELL study (ACTRN12617001579381) who received a donor-derived 4-1BB CAR product targeting CD19. Following infusion, CAR+ T cells and CAR- T cells shows similar differentiation profiles with clonally expanded populations across heterogeneous phenotypes, demonstrating clonal lineages and phenotypic plasticity. We validate these findings in 31 patients with large B cell lymphoma treated with CD19 CAR T therapy. For these patients, we identify using longitudinal mass-cytometry data an association between NK-like subsets and clinical outcomes at 6 months with both CAR+ and CAR- T cells. These results suggest that non-CAR-derived signals can provide information about patients' immune recovery and be used as correlate of clinically relevant parameters.

Mechanisms of mRNA processing defects in inherited THOC6 intellectual disability syndrome.

THOC6 is the genetic basis of autosomal recessive THOC6 Intellectual Disability Syndrome (TIDS). THOC6 facilitates the formation of the Transcription Export complex (TREX) tetramer, composed of four THO monomers. The TREX tetramer supports mammalian mRNA processing that is distinct from yeast TREX dimer functions. Human and mouse TIDS model systems allow novel THOC6-dependent TREX tetramer functions to be investigated. Biallelic loss-of-functon(LOF) THOC6 variants do not influence the expression and localization of TREX members in human cells, but our data suggests reduced binding affinity of ALYREF. Impairment of TREX nuclear export functions were not detected in cells with biallelic THOC6 LOF. Instead, mRNA mis-splicing was observed in human and mouse neural tissue, revealing novel insights into THOC6-mediated TREX coordination of mRNA processing. We demonstrate that THOC6 is required for regulation of key signaling pathways in human corticogenesis that dictate the transition from proliferative to neurogenic divisions that may inform TIDS neuropathology.

Certain heterozygous variants in the kinase domain of the serine/threonine kinase NEK8 can cause an autosomal dominant form of polycystic kidney disease.

Autosomal dominant polycystic kidney disease (ADPKD) resulting from pathogenic variants in PKD1 and PKD2 is the most common form of PKD, but other genetic causes tied to primary cilia function have been identified. Biallelic pathogenic variants in the serine/threonine kinase NEK8 cause a syndromic ciliopathy with extra-kidney manifestations. Here we identify NEK8 as a disease gene for ADPKD in 12 families. Clinical evaluation was combined with functional studies using fibroblasts and tubuloids from affected individuals. Nek8 knockout mouse kidney epithelial (IMCD3) cells transfected with wild type or variant NEK8 were further used to study ciliogenesis, ciliary trafficking, kinase function, and DNA damage responses. Twenty-one affected monoallelic individuals uniformly exhibited cystic kidney disease (mostly neonatal) without consistent extra-kidney manifestations. Recurrent de novo mutations of the NEK8 missense variant p.Arg45Trp, including mosaicism, were seen in ten families. Missense variants elsewhere within the kinase domain (p.Ile150Met and p.Lys157Gln) were also identified. Functional studies demonstrated normal localization of the NEK8 protein to the proximal cilium and no consistent cilia formation defects in patient-derived cells. NEK8-wild type protein and all variant forms of the protein expressed in Nek8 knockout IMCD3 cells were localized to cilia and supported ciliogenesis. However, Nek8 knockout IMCD3 cells expressing NEK8-p.Arg45Trp and NEK8-p.Lys157Gln showed significantly decreased polycystin-2 but normal ANKS6 localization in cilia. Moreover, p.Arg45Trp NEK8 exhibited reduced kinase activity in vitro. In patient derived tubuloids and IMCD3 cells expressing NEK8-p.Arg45Trp, DNA damage signaling was increased compared to healthy passage-matched controls. Thus, we propose a dominant-negative effect for specific heterozygous missense variants in the NEK8 kinase domain as a new cause of PKD.

SRSF1 haploinsufficiency is responsible for a syndromic developmental disorder associated with intellectual disability.

SRSF1 (also known as ASF/SF2) is a non-small nuclear ribonucleoprotein (non-snRNP) that belongs to the arginine/serine (R/S) domain family. It recognizes and binds to mRNA, regulating both constitutive and alternative splicing. The complete loss of this proto-oncogene in mice is embryonically lethal. Through international data sharing, we identified 17 individuals (10 females and 7 males) with a neurodevelopmental disorder (NDD) with heterozygous germline SRSF1 variants, mostly de novo, including three frameshift variants, three nonsense variants, seven missense variants, and two microdeletions within region 17q22 encompassing SRSF1. Only in one family, the de novo origin could not be established. All individuals featured a recurrent phenotype including developmental delay and intellectual disability (DD/ID), hypotonia, neurobehavioral problems, with variable skeletal (66.7%) and cardiac (46%) anomalies. To investigate the functional consequences of SRSF1 variants, we performed in silico structural modeling, developed an in vivo splicing assay in Drosophila, and carried out episignature analysis in blood-derived DNA from affected individuals. We found that all loss-of-function and 5 out of 7 missense variants were pathogenic, leading to a loss of SRSF1 splicing activity in Drosophila, correlating with a detectable and specific DNA methylation episignature. In addition, our orthogonal in silico, in vivo, and epigenetics analyses enabled the separation of clearly pathogenic missense variants from those with uncertain significance. Overall, these results indicated that haploinsufficiency of SRSF1 is responsible for a syndromic NDD with ID due to a partial loss of SRSF1-mediated splicing activity.

Al-Gazali Skeletal Dysplasia Constitutes the Lethal End of ADAMTSL2-Related Disorders.

Lethal short-limb skeletal dysplasia Al-Gazali type (OMIM %601356), also called dysplastic cortical hyperostosis, Al-Gazali type, is an ultra-rare disorder previously reported in only three unrelated individuals. The genetic etiology for Al-Gazali skeletal dysplasia has up until now been unknown. Through international collaborative efforts involving seven clinical centers worldwide, a cohort of nine patients with clinical and radiographic features consistent with short-limb skeletal dysplasia Al-Gazali type was collected. The affected individuals presented with moderate intrauterine growth restriction, relative macrocephaly, hypertrichosis, large anterior fontanelle, short neck, short and stiff limbs with small hands and feet, severe brachydactyly, and generalized bone sclerosis with mild platyspondyly. Biallelic disease-causing variants in ADAMTSL2 were detected using massively parallel sequencing (MPS) and Sanger sequencing techniques. Six individuals were compound heterozygous and one individual was homozygous for pathogenic variants in ADAMTSL2. In one of the families, pathogenic variants were detected in parental samples only. Overall, this study sheds light on the genetic cause of Al-Gazali skeletal dysplasia and identifies it as a semi-lethal part of the spectrum of ADAMTSL2-related disorders. Furthermore, we highlight the importance of meticulous analysis of the pseudogene region of ADAMTSL2 where disease-causing variants might be located. © 2023 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

DNA methylation episignature for Witteveen-Kolk syndrome due to SIN3A haploinsufficiency.

PURPOSE: Witteveen-Kolk syndrome (WITKOS) is a rare, autosomal dominant neurodevelopmental disorder caused by heterozygous loss-of-function alterations in the SIN3A gene. WITKOS has variable expressivity that commonly overlaps with other neurodevelopmental disorders. In this study, we characterized a distinct DNA methylation epigenetic signature (episignature) distinguishing WITKOS from unaffected individuals as well as individuals with other neurodevelopmental disorders with episignatures and described 9 previously unpublished individuals with SIN3A haploinsufficiency. METHODS: We studied the phenotypic characteristics and the genome-wide DNA methylation in the peripheral blood samples of 20 individuals with heterozygous alterations in SIN3A. A total of 14 samples were used for the identification of the episignature and building of a predictive diagnostic biomarker, whereas the diagnostic model was used to investigate the methylation pattern of the remaining 6 samples. RESULTS: A predominantly hypomethylated DNA methylation profile specific to WITKOS was identified, and the classifier model was able to diagnose a previously unresolved test case. The episignature was sensitive enough to detect individuals with varying degrees of phenotypic severity carrying SIN3A haploinsufficient variants. CONCLUSION: We identified a novel, robust episignature in WITKOS due to SIN3A haploinsufficiency. This episignature has the potential to aid identification and diagnosis of individuals with WITKOS.

X-Linked intellectual disability update 2022.

Genes that are involved in the transcription process, mitochondrial function, glycoprotein metabolism, and ubiquitination dominate the list of 21 new genes associated with X-linked intellectual disability since the last update in 2017. The new genes were identified by sequencing of candidate genes (2), the entire X-chromosome (2), the whole exome (15), or the whole genome (2). With these additions, 42 (21%) of the 199 named XLID syndromes and 27 (25%) of the 108 numbered nonsyndromic XLID families remain to be resolved at the molecular level. Although the pace of discovery of new XLID genes has slowed during the past 5 years, the density of genes on the X chromosome that cause intellectual disability still appears to be twice the density of intellectual disability genes on the autosomes.

High activation levels maintained in receptor-binding domain-specific memory B cells in people with severe coronavirus disease 2019.

The long-term health consequences of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection are still being understood. The molecular and phenotypic properties of SARS-CoV-2 antigen-specific T cells suggest a dysfunctional profile that persists in convalescence in those who were severely ill. By contrast, the antigen-specific memory B-cell (MBC) population has not yet been analyzed to the same degree, but phenotypic analysis suggests differences following recovery from mild or severe coronavirus disease 2019 (COVID-19). Here, we performed single-cell molecular analysis of the SARS-CoV-2 receptor-binding domain (RBD)-specific MBC population in three patients after severe COVID-19 and four patients after mild/moderate COVID-19. We analyzed the transcriptomic and B-cell receptor repertoire profiles at ~2 months and ~4 months after symptom onset. Transcriptomic analysis revealed a higher level of tumor necrosis factor-alpha (TNF-α) signaling via nuclear factor-kappa B in the severe group, involving CD80, FOS, CD83 and TNFAIP3 genes that was maintained over time. We demonstrated the presence of two distinct activated MBCs subsets based on expression of CD80hi TNFAIP3hi and CD11chi CD95hi at the transcriptome level. Both groups revealed an increase in somatic hypermutation over time, indicating progressive evolution of humoral memory. This study revealed distinct molecular signatures of long-term RBD-specific MBCs in convalescence, indicating that the longevity of these cells may differ depending on acute COVID-19 severity.

Clinical and functional characterization of germline PIK3CA variants in patients with PIK3CA-related overgrowth spectrum disorders.

Mosaic variants in the PIK3CA gene, encoding the catalytic subunit of phosphoinositide 3-kinase (PI3K), produce constitutive PI3K activation, which causes PIK3CA-related overgrowth spectrum disorders. To date, fewer than 20 patients have been described with germline alterations in PIK3CA. In this study, we describe three unrelated individuals with overgrowth and germline PIK3CA variants. These variants were discovered through whole-exome sequencing and confirmed as germline by testing multiple tissue types, when available. Functional analysis using Patient 1's fibroblast cell line and two previously reported patients' cell lines showed increased phosphorylation of AKT during cellular starvation revealing constitutive activation of the phosphoinositide-3-kinase/protein kinase B/mechanistic target of rapamycin (PI3K/AKT/mTOR) pathway. Alternatively, stimulation of the cells by fetal bovine serum produced a reduced response, indicating an activated status of the PI3K complex reducing the pathway response to further external stimulation. Additional studies utilizing Biolog Phenotype Microarray technology indicated reduced energy production when cells were exposed to growth factors stimulating the PI3K/AKT/mTOR pathway, confirming the trend observed in the AKT phosphorylation test after stimulation. Furthermore, treatment with inhibitors of the PI3K/AKT/mTOR pathway rescued the normal energy response in the patients' cells. Collectively, these data demonstrate that disease-causing germline PIK3CA variants have a functional consequence, similar to mosaic variants in the PI3K/AKT/mTOR pathway.

Identification of human progenitors of exhausted CD8+ T cells associated with elevated IFN-γ response in early phase of viral infection.

T cell exhaustion is a hallmark of hepatitis C virus (HCV) infection and limits protective immunity in chronic viral infections and cancer. Limited knowledge exists of the initial viral and immune dynamics that characterise exhaustion in humans. We studied longitudinal blood samples from a unique cohort of individuals with primary infection using single-cell multi-omics to identify the functions and phenotypes of HCV-specific CD8+ T cells. Early elevated IFN-γ response against the transmitted virus is associated with the rate of immune escape, larger clonal expansion, and early onset of exhaustion. Irrespective of disease outcome, we find heterogeneous subsets of progenitors of exhaustion, based on the level of PD-1 expression and loss of AP-1 transcription factors. Intra-clonal analysis shows distinct trajectories with multiple fates and evolutionary plasticity of precursor cells. These findings challenge the current paradigm on the contribution of CD8+ T cells to HCV disease outcome and provide data for future studies on T cell differentiation in human infections.

Inferring Epistasis from Genetic Time-series Data.

Epistasis refers to fitness or functional effects of mutations that depend on the sequence background in which these mutations arise. Epistasis is prevalent in nature, including populations of viruses, bacteria, and cancers, and can contribute to the evolution of drug resistance and immune escape. However, it is difficult to directly estimate epistatic effects from sampled observations of a population. At present, there are very few methods that can disentangle the effects of selection (including epistasis), mutation, recombination, genetic drift, and genetic linkage in evolving populations. Here we develop a method to infer epistasis, along with the fitness effects of individual mutations, from observed evolutionary histories. Simulations show that we can accurately infer pairwise epistatic interactions provided that there is sufficient genetic diversity in the data. Our method also allows us to identify which fitness parameters can be reliably inferred from a particular data set and which ones are unidentifiable. Our approach therefore allows for the inference of more complex models of selection from time-series genetic data, while also quantifying uncertainty in the inferred parameters.

Mosaicism of common pathogenic MECP2 variants identified in two males with a clinical diagnosis of Rett syndrome.

Rett (RTT) syndrome, a neurodevelopmental disorder caused by pathogenic variation in the MECP2 gene, is characterized by developmental regression, loss of purposeful hand movements, stereotypic hand movements, abnormal gait, and loss of spoken language. Due to the X-linked inheritance pattern, RTT is typically limited to females. Recent studies revealed somatic mosaicism in MECP2 in male patients with RTT-like phenotypes. While detecting mosaic variation using Sanger sequencing is theoretically possible for mosaicism over ~15%-20%, several variables, including efficiency of PCR, background noise, and/or human error, contribute to a low detection rate using this technology. Mosaic variants in two males were detected by next generation sequencing (NGS; Case 1) and by Sanger re-sequencing (Case 2). Both had targeted digital PCR (dPCR) to confirm the variants. In this report, we present two males with classic RTT syndrome in whom we identified pathogenic variation in the MECP2 gene in the mosaic state (c.730C > T (p.Gln244*) in Patient 1 and c.397C > T (p.Arg133Cys) in Patient 2). In addition, estimates and measures of mosaic variant fraction were surprisingly similar between Sanger sequencing, NGS, and dPCR. The mosaic state of these variants contributed to a lengthy diagnostic odyssey for these patients. While NGS and even Sanger sequencing may be viable methods of detecting mosaic variation in DNA or RNA samples, applying targeted dPCR to supplement these sequencing technologies would provide confirmation of somatic mosaicism and mosaic fraction.